Pair supersolid of the extended Bose-Hubbard model with atom-pair hopping on the 

triangular Lattice 
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We systematically study an extended Bose-Hubbard model with atom hopping and atom-pair 
hopping in the presence of a three-body constraint on the triangular lattice. By means of large- 
scale Quantum Monte Carlo simulations, the ground-state phase diagram are studied. We find a 
continuous transition between the atomic superfluid phase and the pair superfluid when the ratio of 
the atomic hopping and the atom-pair hopping is adapted. We then focus on the interplay among the 
atom-pair hopping, the on-site repulsion and the nearest-neighbor repulsion. With on-site repulsion 
present, we observe first order transitions between the Mott Insulators and pair superfluid driven by 
the pair hopping. With the nearest-neighbor repulsion turning on, three typical solid phases with 
2/3, 1 and 4/3-filling emerge at small atom-pair hopping region. A stable pair supersolid phase is 
found at small on-site repulsion. This is due to the three-body constraint and the pair hopping, which 
essentially make the model a quasi hardcore boson system. Thus the pair supersolid state emerges 
basing on the order-by-disorder mechanism, by which hardcore bosons avoid classical frustration on 
the triangular lattice. The transition between the pair supersolid and the pair superfluid is first 
order, except for the particle- hole symmetric point. We compare the results with those obtained by 
means of mean-field analysis. 

PACS numbers: 75.10.Jm, 05.30. Jp, 03.75.Lm, 37.10.Jk 



I. INTRODUCTION 

The development in experimentally manipulating 
ultra-cold atoms in an optical lattice paved the way 
to simulate strongly interacting systems in condensed- 
matter physics^. It provides very clean and tunable 
systems to study quantum phase transitions and exotic 
quantum states, e.g., superfluid, Mott insulator, which 
are not easily accessible in condensed matters. The con- 
densation of paired electrons, which provides the basis of 
superconductivity of metallic superconductor, plays an 
essential role in modern condensed-matter physics. Thus 
realizing pairing related novel quantum states in the con- 
text of ultra-cold atoms has attracted considerable recent 
interest, both in theoretical and experimental research. 

One candidate to achieve such states is lattice bosons 
with attractive on-site interactions, which is stabilized by 
a three-body constraint^. The three-body constraint 
has been realized by large three-body loss processes^. 
The system can be mapped into spin-one atoms at unit 
filling*. Besides the conventional single atom superfluid 
(ASF) state, a pair (dimer) superfluid (PSF) phase con- 
sisting of the condensation of boson pairs emerges under 
sufficiently strong attraction 2 -^. The PSF state is man- 
ifested as a second-order effect in the atom hopping in 
the optical lattice. Various phase transitions among the 
atomic superfluid (ASF), Mott insulator (MI) and PSF 
are investigated in great detail 2 -^ - — . Both the ground 
state and the thermal phase diagrams are obtained. Al- 
though the pair supersolid (PSS) state was predicted in 
system with correlated hoppin g 12 ' 13 and the paired two- 
species bosons supersolid was predicted on the square 



lattice^ with the attractive interaction between differ- 
ent species and the repulsion between the same species 
of atoms turning on, the pair (single species) supersolid 
was not found in the present system when the nearest- 
neighbor repulsion is included, except for an isolated con- 
tinuous supersolid at the Dirac point 2 -. The reason might 
be the same instability of the supersolid state (SS) on the 
square lattice, on which the former research focused, for 
hardcore boson a 15 ' 16 . 

Another practical way to access pairing phenomena is 
considering repulsive bound-atom pairs in optical lattice, 
which have been realized in experiments^^. The system 
can be reasonably described by explicitly including atom- 
pair hoppin g 19 ' 20 in the ordinary extended Bose-Hubbard 
model (EBH). Such a pair hopping can also be introduced 
in atom-molecule coupling system on the state-dependent 
optical lattice^, or by a mechanism based on transport- 
inducing collisions^ 2 -. The resulted Bose-Hubbard Hamil- 
tonian is thus 
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h = - (Mi + — n i( n i ~ !) + Vn 
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where aj (a*) creates (annihilates) a boson at site i, t is 
the atom hopping amplitude for nearest neighbor sites 
(ij), t p is the pair-hopping amplitude, V the nearest- 
neighbor repulsion, fi the chemical potential, U > the 
on-site repulsion. The three-body constraint is applied 
which requires the occupation number m — 0,1, or 2. 
It is naturally expected that PSF, MI and solid phases 
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should emerge at certain parameters. 

The three-body constraint and the pair hopping makes 
the system resemble hardcore bosons. For hardcore 
bosons supersolid state emerges on the triangular lattice 
basing on an order-by-disorder mechanism, by which a 
quantum system avoids classical frustratio n 23 ^ 4 . Aim- 
ing to realize PSS state, we thus focus on the triangular 
lattice. A preliminary mean-field (MF) analysis 2 ^ does 
predict the PSS phase in such a system. In present work 
we systematically study the model by means of a large 
scale Quantum Monte Carlo simulation. The phase dia- 
gram of the system is studied in great detail. With the 
nearest-neighbor interaction turning on, we report three 
types of solid phases at small hopping strength. Increas- 
ing pair hopping t p , we find the expected PSS phase. 

This paper is organized as follows: We first discuss the 
ground state at the classical limit in Sec[IIl which is use- 
ful to identify various solids and MI states, in which we 
are interested. We describe the quantum Monte Carlo 
(QMC) method and useful observables in Sec. IIIII We 
then present, in Sec lIVl QMC simulation results, com- 
paring with the mean- field analysis. The results for the 
noninteracting (U = V = 0) case are described in Sec. 
IIV A[ Those for the interacting case U ^ 0, V = are 
given in Sec. IIVBI and for the case V ^ in Sec. IIV CI 
focusing on the parameter regions where the pair super- 
solid phase emerges. The the nature of related phase 
transitions are discussed. We conclude in Sec. |V] with 
discussions. 



II. CLASSICAL LIMIT 




FIG. 1. (Color online) (a) The ground-state (U/V - 
fi/V) phase diagram in the classical limit t = 0, t p = 0. 
(nA,riB, nc) denotes a solid order in a y3 x y3 ordering with 
wave vector Q = (4-7r/3, 0), where the numbers are boson oc- 
cupations on three sublattices, respectively, (b) SI phase (0, 
0, 2) with filling p = 2/3. (c) SII phase (0, 1, 2) with filling 
p = 1. (d) SIII phase (0, 2, 2) with filling p = 4/3. 



We first study the zero temperature phase diagram in 
the classical limit (t = 0, t p = 0) in the presence of the 
three-body constraint. 

At half filling, the solid order is frustrated on the trian- 
gular lattice, and the classical model has a hugely degen- 
erate ground state with an extensive zero-temperature 
entropy 2 ^. We focus on the solid ordering, away from 
half filling, with wave vector Q = (47r/3,0). The lattice 
is divided into three sublattices: A, B, C. The solid order 
is thus denoted as (jiA,n B ,nc) or its equivalent permu- 
tations, where the numbers are the boson occupations on 
three sublattices, respectively. The energy per site is 

e = -fi h V{n A n B + n A n c + n B n c ) 

+ jj iT-AjriA - 1) + n B (n B - 1) + n c (n c - 1) ,^ 
6 

By comparing energy per site, we obtain phase bound- 
aries between various phases, as shown in Fig. [1] Three 
atom-pair solid phases appearing in the region < 
U/V < 3 are of special interest, which have two bosons 
sitting on one or two sublattices, as shown in FigQJb), 
(c) and (d). They are denoted as SI, SII and SIII, with 
densities p = 2/3, 1,4/3, respectively. Other solid states 
related with configurations (0, 0, 1), (1, 1, 0), (1, I, 2), 



(1, 2, 2) and MI states (1, 1, 1) and (2, 2, 2) are also 
found. 



III. METHOD AND OBSERVABLES 

We simulate the model JT]) using the stochastic series 
expansion (SSE) QMC method^! with the directed loop 
update 28 , which is improved for the pair hopping. The 
head of a directed loop carries a create (annihilate) op- 
erator a(a^) in the conventional directed loop algorithm 
for the BH model with single particle hopping. We im- 
prove the algorithm by allowing the head of a directed 
loop carry a pair create (annihilate) operators a 2 ' (a 2 ) 
for the present EBH model with additional pair hopping 
terms. Similar improvement has been described in the 
literature 2 ^. In the simulations, temperature is chosen as 
(3 = L to ensure the system sitting in its ground state. 

To distinguish the ASF and the PSF states, we de- 
fine an even (odd) superfluid stiffness as the order 
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where (3 is the inverse temperature, a can be 'even' or 
'odd'; W(a) is the total even (odd) winding number. The 
factor '4' multiplying t p is due to the hopping of paired 

atoms^i. The total superfluid stiffness is p s = pi even ^ + 
pi° dd \ For a PSF state, we define the pair superfluid 
order parameter pf^ = pf ven ^ — p < f dd * > > 0, while for an 
ASF state, p^ = 0, but p s ^ 0. It is worthy to note 
that, without single atom hopping, p~° dd ^ is guaranteed 
to be zero. 

The structure factor is defined to characterize the solid 
order: 

S(Q)/JV=(pqpJ,>> (4) 

where pq — 0-/N) exp(iQri). The solid states SI, 

SII and SIII share the same ordering at the wave vector 
Q = (4tt/3,0), but bear different values 4/9, 1/3, 4/9, 
respectively, in the perfect ordering. 

IV. RESULTS 

A. Non-interaction case: U = 0,V = 0; Competition 
between ASF and PSF 



particle picture: to put a boson on the empty lattice 
gains potential energy — p, and kinetic energy — zt, where 
z = 6 is the coordination number of the triangular lattice. 
Adding a hole to the MI (p = 2) state costs chemical en- 
ergy pi, while obtain kinetic energy — 2zt. At large pair 
hopping t p ^S> t, the system shows a PSF-Empty tran- 
sition at fi/t = —Qtp/t. A pair of bosons emerging on 
the Empty phase gets —2zt p kinetic energy and poten- 
tial energy —2//. Similarly, the PSF-MI transition line 
is at ji/t = 6t p /t by analyzing the emerging of a pair of 
holes on the MI state. 

To demonstrate the transition between the ASF and 
the PSF states in more detail, we show here the system 
behaviors along p = in Fig. Hfb). The PSF state was 
predicted at t p /t > 2.9 in the MF frame using the crite- 
rion (j>i = (a) = but 02 = (a 2 ) ^ 02£. This transition is 
confirmed by our unbiased QMC simulations. However, 
the transition point is t p /t — 2.5. The PSF state is char- 
acterized by pi p ^ > 0. The phase transition between ASF 
and PSF is continuous. In the presence of strong on-site 
attractive interactions (U < 0), the PSF has been pre- 
dicted in the lattice bosons with a three-body hard-core 
constraint 8 ' 9 . However, the transition between ASF and 
PSF phases is claimed to be first order—. 




B. MI-PSF phase transitions for U > and V = 
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FIG. 2. (Color online) (a) (zt p /t, p/t) phase diagram of model 
(TTJ) obtained by means of the MF (solid blue lines) and the 
QMC(circles) method, which contains Empty, MI, SF, and 
PSF phases at U = 0, V = 0. (b) Pair superfluid order 
parameter ps as a function of zt p /t for p = (L = 12). The 
MF order parameters <f>i and <j>2 are also shown. 



The ASF-PSF transition for model (JT]) without on-site 
repulsion U and nearest-neighbor repulsion V was pre- 
viously discussed by means of MF analysis^. Here we 
present our QMC results for such a system. To demon- 
strate the differences between our results and the MF 
results, we show both of them in Fig. [2J We find an 
Empty phase and a p = 2 MI phase at negative and large 
chemical potentials, respectively. The presence of the MI 
state is due to the three-body constraint. Between the 
Empty phase and the MI phase, there are two SF phases: 
an ASF phase and a PSF phase. 

With small pair hopping t p ~ 0, the system exhibits 
an Empty- ASF transition at p = — 6t and a SF-MI tran- 
sition at p — 12t. This can be understood in the single 
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FIG. 3. (Color online) (zt p /U,p/U) phase diagram ob- 
tained by QMC (circles) and MF analysis (solid blue lines) 
at t = 0, V = 0, which contains the Empty, MI (p = 1), 
MI (p = 2), and PSF phases at p/V = 0.55. Open and 
solid circles correspond to continuous and first order transi- 
tions, respectively. The pair superfluid order parameter p^ 
vs. ztp/U is shown in the inset. 

In the presence of the on-site repulsion U and the ab- 
sence of the single atom hopping t, the ASF state dis- 
appears, instead an MI(p = 1) phase emerges. Without 
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three-body constraint, the general picture of the phase 
diagram has been studied by Zhou et al using the MF 
method^. We show here the phase diagram, under three- 
body constraint, obtained by using the MF method— and 
by the QMC simulations in Fig. [3l The phase boundaries 
obtained by using the two methods are in good agree- 
ment for large t P . The Empty-PSF phase transition line 
is straight with the slope p/(zt p ) = — 1, which can be 
understood in the single particle picture: A boson pair 
emerging in the empty phase gets —2zt p kinetics energy 
and — 2/i potential energy. Similarly, the MI(p = 2)- 
PSF boundary is also straight with a slope p/(zt p ) = 1 
at large t p . QMC simulations provide more accurate 
boundaries between the MI phases and PSF phase. The 
p = 2 MI-PSF transition is found to be continuous and 
the MI(p = 1)-PSF transition is first order, as illustrated 
in the inset of Fig. [3] 



and QMC simulations is shown in Fig [4] Both method 
predict a PSS phase between the two solid phases. How- 
ever, the PSS region is much smaller in the QMC phase 
diagram. This is due to quantum fluctuations ignored in 
the MF analysis. 

We further illustrate the phase diagram along two 
lines: one is zt p /V = 1.05, the other is p/V = 7.5. 

Figure [5] shows the density p, structure factor S(Q)/N 

(p) 

and pair superfluid order parameter p s as functions of 
p/V along the zt p /V = 1.05 line. Considering the hole- 
particle symmetry, we only present results for p/V < 
6. The PSS state emerges clearly between the two solid 
phases (SI and SIII), which is further proved by a finite- 
size scaling analysis of the structural factor and the pair 

superfluid stiffness pi (not shown). 



Pair supersolid phase for U — and U > 




FIG. 4. (Color online) (zt p /V, p/V) phase diagram obtained 
by QMC (circles) and the MF analysis (solid blue lines) of 
model (ff]) at t — 0, U = 0, which contains the Empty, MI, SI, 
SIII, PSF and PSS phases. Open and solid circles correspond 
to continuous and first order transitions, respectively. The 
PSS-PSF transition is first order except for the p/V = 6 
point marked by a red square^. 

As a result of the absence of on-site repulsion, boson 
pairs easily emerge on a site without costing of any po- 
tential energy. In the limit t p /t — > oo (or t = 0), boson 
pairs (hole pairs) can easily hop on the lattice, form- 
ing a PSF phase. When the nearest-neighbor interaction 
V present, solid ordering appears naturally. Hardcore 
bosons with single-atom hopping on the triangular lat- 
tice has an SS phase between two solid phases ((0, 0, 1) 
and (0, 1, l)) 23 i 24 . Bosons on the triangular lattice with 
pair hopping under three-body constraint mimic such be- 
havior with the two solid phases replaced by the (0, 0, 
2) and (0, 2, 2) solid phases respectively. The global 
phase diagram obtained by means of the MF method 25 
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FIG. 5. (Color online) The density p, structural factor 
S(Q)/N, PSF stiffness pi p) as functions of p/V at zt p = 1.05. 

Next, we scan the pair hopping t p along p/V — 7.5. 
The density p, structural factor S(Q)/N and pair super- 
fluid order parameter ps as functions of the pair hopping 
t p are shown in Fig. EJa), with the PSS state emerging 
in the region 1.668 < zt p /V < 1.788. Figure EJb) shows 
a finite size scaling analysis of the pair superfluid or- 
der parameter and the structural factor S(Q)/N at 
ztp/V — 1.74, p/V = 7.5, which proves that the state is 
indeed a pair supersolid. 

The phase transition between the PSS and other states 
are interesting. The PSS-SI(SIII) transition is continuous 
from continuous variation of the density and the struc- 
ture factor, while the PSS-PSF transition is first order, 
except for the particle-hole symmetric points, as illus- 
trated by the jump of the structural factor and the PSF 
order parameter in Fig. [5{a). From the jumps of the PSF 
order parameter and the structural factor (not shown), 
the PSF-solid phase transitions are also first order. 

We now check whether the pair supersolid phase could 
survive as a weak on-site repulsion turns on. Figure [7] 
shows the MF and QMC phase diagram at weak repul- 
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FIG. 6. (Color online) (a) The density p, structural factor 
S(Q)/N, and PSF order parameter p^ as functions of pair 
hopping ztp/V at fi/V = 7.5. (b) Finite size scalings of the 
structural factor S(Q)/N and PSF order parameter pf 1 at 
ztp/V — 1.74 and p/V — 7.5, showing PSS character at ther- 
modynamic limit. 



sion U/V = 1. Three new solid phases, i.e., SII, (0, 0, 1) 
and (2, 2, 1), emerge comparing to the case U = 0. The 
PSS phase persists. Increasing t p , all solid phases melt 
into the PSF phase eventually. The persistence of the 
PSS state at weak on-site repulsion is because that the 
PSS state depends strongly on the frustrated geometri- 
cal structure of the triangular lattice for hardcore bosons: 
although the maximum occupation is two, the pair hop- 
ping essentially make the model ([1} a quasi hardcore bo- 
son system. Weak on-site repulsions do not change this 
feature. 

At the large U/V limit, the phase diagram should re- 
duce to the one shown in Fig. [31 This suggests the PSS 
phase will disappear finally. Indeed, we don't find PSS 
phase when U/V > 4. 




FIG. 7. (Color online) (zt p /V, p/V) phase diagram obtained 
by QMC (circles) and the MF (solid blue lines) method at 
U/V = 1, which contains the Empty, Ml(p = l,p = 2), SI, 
SII, SIII, (0, 0, 1) , (2, 2, 1), PSF and PSS phases. Open and 
solid circles correspond to continuous and first order transi- 
tions, respectively. The red square marks the special high 
symmetric point. 



V. CONCLUSION 

We have studied the ground-state phase diagram of the 
extended Bose-Hubbard model with explicit atom-pair 
hopping terms, under the three-body constraint, on the 
triangular lattice. Experimentally, the atom-pair hop- 
ping physics has been realized. The triangular optical lat- 
tice can also be implemented by three beams of lasers 33 . 
The atom-molecule coupling system^! and the spin-1 bo- 
son system^ in state dependent optical triangular lattice 
are also candidates to simulate pair related quantum phe- 
nomena. By means of the improved SSE QMC method, 
we obtain accurate phase diagrams at various parameter 
values. Rich solid phases and the pair supcrfluid phase 
are found, as well as the pair supersolid phase, which 
emerges when the nearest-neighbor repulsion and the pair 
hopping are present, under the condition that the on-site 
repulsion is not large. The three-body constraint and 
the pair hopping essentially make the model ([T]) a quasi 
hardcore boson system. Weak on-site repulsions do not 
change this feature. Therefore, the mechanism of form- 
ing the PSS state is the same as that of forming SS state 
for hardcore bosons on the triangular lattice. The prop- 
erties of phase transitions involved, e.g., the ASF-PSF, 
the MI-PSF, the PSS-PSF, the PSF-solid transitions, are 
studied. 

It would be interesting to study the finite temperature 
properties of the present model. The three-body con- 
straint due to the three-body loss^, can be suppressed, 
or even inhibited by the quantum Zeno effecli^. Thus 
the quantum Zeno effect provides us a basis to study the 
pair physics without three-body constraint. 
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